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Abstract: We extend to rank-based tests of multivariate independence the 
Chcrnoff-Savage and Hodges-Lehmann classical univariate results. More pre- 
cisely, we show that the Taskinen, Kankainen and Oja (2004) normal-score rank 
test for multivariate independence uniformly dominates -in the Pitman sense — 
the classical Wilks (1935) test, which establishes the Pitman non-admissibility 
of the latter, and provide, for any fixed space dimensions p, q of the marginals, 
the lower bound for the asymptotic relative efficiency, still with respect to 
Wilks' test, of the Wilcoxon version of the same. 

1. Introduction 

1.1. Testing multivariate independence 

The problem of testing for independence between two random variables with un- 
specified densities has been among the very first applications of ranks in statistical 
inference. Spearman's correlation coefficient was proposed as early as 1904 (Spear- 
man [29]), long before Wilcoxon [33] 's rank sum and signed rank tests for location, 
and Kendall introduced his rank correlation measure in 1938 (Kendall [18]). 

The multivariate version of the same problem, namely, testing independence be- 
tween two random vectors with unspecified densities, is harder. The first rank-based 
counterpart to the Gaussian likelihood ratio method of Wilks [34] was developed 
in Chapter 8 of the monograph by Puri and Sen [25] and, for almost thirty years, 
has remained the only rank-based solution to the problem. The proposed test, how- 
ever, is based on componentwise rankings, hence is neither invariant under affinc 
transformations, nor distribution- free - unless of course both random vectors have 
dimension one, in which case we are back to the traditional context of rank-based 
tests of bivariate independence. 

Of course, the asymptotic null distributions of the Puri and Sen test statistics 
do not depend on the underlying distributions. This however is a simple conse- 
quence of central-limit behavior and studentization - the matrix of the quadratic 
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form involved being a consistent empirical version of its population counterpart. 
Such property, in the terminology of Randies [27] is only a weak form of asymp- 
totic distribution-freeness, the strong form being asymptotic equivalence with a 
genuinely distribution-free random variable. 

Alternatives to the Puri and Sen test have been made possible by the recent 
developments of multivariate concepts of signs and ranks. Gieser and Randies [3] 
propose a sign test based on Randies [28] 's inter directions. Taskinen, Kankainen and 
Oja [30] also propose a sign test, based on the so-called standardized spatial signs, 
which is asymptotically equivalent to the Gieser and Randies one under elliptic 
symmetry assumptions. Multivariate ranks are introduced, along with the signs, in 
Taskinen, Oja and Randies [32], where multivariate analogs of Spearman's rho and 
Kendall's tau (based on Wilcoxon scores) are considered. This is extended, in Task- 
inen, Kankainen and Oja [31], to arbitrary score functions. All these tests rely on 
exactly affinc- invariant and strongly (in Randies' sense) asymptotically distribution- 
free statistics - this latter property, however, is obtained at the expense of stronger 
assumptions of elliptical symmetry. 

The objective of this paper is to show that Chernoff-Savage and Hodges-Lehmann 
results can be established for the normal (van der Waerden, say) and Wilcoxon score 
versions of the Taskinen, Kankainen and Oja [31] tests, establishing their excellent 
performances with respect to the Gaussian likelihood ratio procedure of Wilks [34] . 

1.2. Chernoff-Savage and Hodges-Lehmann results 

Rank-based inference long has been considered as a somewhat heuristic and het- 
eroclite collection of "quick-and-easy" methods applicable under a broad range of 
assumptions. The Spearman, Kendall, Wilcoxon, . . . tests were not introduced in 
the context of a general theory, but as isolated tools, to be used when everything 
else fails. Certainly, their performances were not expected to compare favorably 
with those of their parametric competitors. 

This widespread opinion was partly dispelled by two famous papers - Hodges and 
Lehmann [15] and Chcrnoff and Savage [1] -establishing that rank-based methods 
not only compete very well, but even may outperform their parametric counter- 
parts (quite remarkably, Hotelling and Pabst [17] already show that the asymptotic 
relative efficiency, under bivariate Gaussian densities, of Spearman with respect to 
the traditional correlation coefficient, is as high as 0.9119). These papers certainly 
played a triggering role in the subsequent development of the well-structured theory 
of rank-based inference associated, mainly, with the name of Hajek, culminating in 
the monographs by Hajek and Sidak (sec Hajek, Sidak and Sen [4]) and Puri and 
Sen [26]. 

In their celebrated ".864 theorem", Hodges and Lehmann [15] proved that, in 
the two-sample location model (but this extends to more general location problems, 
such as one-sample, c-sample, AN OVA, regression, etc.), the Pitman asymptotic 
relative efficiency (ARE) of Wilcoxon (i.e., linear-score) rank tests with respect to 
their normal-theory competitor (namely, Student's standard two-sample i-test) is 
never less than .864. In other words, irrespective of the underlying distribution, 
Wilcoxon tests asymptotically never need more than 13.6% observations more than 
i-tests to achieve the same power. 

No less surprising is the Chernoff and Savage [1] result establishing the amazing 
fact that, in the same class of models, the ARE of van der Waerden (i.e., normal- 
score) rank tests, still with respect to the corresponding standard Gaussian ones, 
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are always larger than one, and that this minimal value is reached at Gaussian 
distributions only. It should be stressed that these results deal with the "worst 
cases" : both for the Wilcoxon and the van der Waerden tests, it is possible to show 
that there is no "best case" , in the sense that it is possible to construct a sequence of 
underlying distributions along which AREs (still with respect to standard Gaussian 
tests) tend to infinity. 

With the extensions of rank-based inference to other models and other inference 
problems, it has been discovered that these Chernoff-Savage and Hodges-Lehmann 
results are not just a happy accident specific to location problems and univariate 
observations. Hallin [5] showed that the van der Waerden version of the serial 
rank tests proposed by Hallin and Puri [13] also uniformly beats (in the Pitman 
sense) the corresponding everyday practice parametric Gaussian tests based on 
autocorrelations. As for the extension of the Hodges-Lehmann [15] result to this 
time series setup, the lower bound of the Wilcoxon version of those tests, still 
with respect to the parametric Gaussian tests, was shown to be .856 by Hallin and 
Tribcl [14]. 

Extensions to (nonserial as well as serial) problems involving multivariate ob- 
servations were recently obtained by Hallin and Paindaveine [7, 8, 10], who showed 
that their various multivariate van der Waerden rank tests uniformly dominate the 
corresponding parametric Gaussian procedures in a broad class of problems (culmi- 
nating in the problem of testing linear restrictions on the parameters of a multivari- 
ate general linear model with vector ARM A errors); the Pitman non-admissibility 
of the associated everyday practice Gaussian tests (one-sample and two-sample 
Hotelling tests, multivariate -F-tests, multivariate portmanteau and Durbin- Watson 
tests, etc.) follows. Hallin and Paindaveine [7] (resp., Hallin and Paindaveine [8]) 
also extended Hodges-Lehmann's result to the multivariate location (resp., serial) 
setup, providing, for any fixed dimension of the observation space, the lower bound 
for the AREs of the proposed multivariate linear-score (i.e., Wilcoxon type) rank 
tests with respect to the parametric Gaussian ones. 

Results recently have also been obtained beyond the case of location parameters 
(in a broad sense, including regression and autoregression coefficients): Paindav- 
eine [23] shows that the uniform Pitman dominance of normal-score rank-based 
procedures over their parametric Gaussian competitors extends to the shape pa- 
rameter of elliptical populations; see Hallin and Paindaveine [11] and Hallin, Oja 
and Paindaveine [6]. It follows that, for any space dimension k > 2, the Gaussian 
maximum likelihood estimator of shape, based on empirical covariances, as well as 
the corresponding Gaussian likelihood ratio tests, are Pitman-nonadmissible. 

In this paper, similar results arc established for the problem of testing inde- 
pendence between elliptical random vectors. More precisely, we consider the van 
der Waerden and Wilcoxon versions of the rank score tests recently proposed by 
Taskinen, Kankainen and Oja [31]. We first establish, under (multivariate exten- 
sions of) the Konijn [19] alternatives, a Chernoff-Savage result, showing that the 
van der Waerden version of their test uniformly dominates the parametric Gaus- 
sian procedure - Wilks [34] 's test - which establishes the Pitman-nonadmissibility of 
the latter. Similarly, we extend the Hodges-Lehmann [15] ".864 result", providing, 
for any fixed couple of marginal space dimensions (p, q) , the lower bound for the 
asymptotic relative efficiency of the Taskinen, Kankainen and Oja [31] Wilcoxon 
test with respect to Wilks. 
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1 . 3. Outline of the paper 

The paper is organized as follows. In Section 2, we define the notation to be used 
throughout, describe the problem of testing for independence between elliptically 
symmetric marginals, and briefly recall the rank score tests developed by Taskinen, 
Kankainen and Oja [31]. In Section 3, we establish the Pitman non-admissibility 
of Wilks' test for multivariate independence. The analog of Hodges-Lehmann [15] 's 
result for the problem under study is derived in Section 4. Section 5 briefly con- 
cludes. 

2. Rank-based tests for multivariate independence 
2.1. Elliptical symmetry 

Recall that the distribution of a random fc-vector X is said to be elliptically sym- 
metric if and only if its probability density function is of the form 

t 2 - 1 ) Z^ S;/ (x):^c fcJ (detS)- 1 /2 / (((x-/x)'E- 1 (x- / x)) 1/2 ), xel 1 , 

for some k- vector /x (the centre of the distribution), some symmetric positive definite 
real k x k matrix S = (S^) with |E| = 1 (the shape matrix), and some function 
/ : M.q — > R + such that / > a.e. and fJ-k-i-j '■= r f(r) dr < oo (ckj is a 
normalization factor depending on the dimension k and /). Denote by Pfe(/x,S, /) 
the corresponding distribution. 

The shape matrix S determines the orientation and shape of the equidensity 
contours associated with / which are a family of hyper-ellipsoids centered at u. 
The definition we are adopting here involves a determinant-based normalization 
which, in a sense, can be considered as canonical : see Paindaveine [24] and Hallin 
and Paindaveine [8] for details. Other normalizations are also possible (such as 
En = 1 or tr£ = fc); the choice of such a normalization does not play any role in 
the sequel. 

The problem of testing for multivariate independence (see Section 2.2 below) 
is invariant under (block-)affine transformations, and so are all tests considered in 
this paper. Therefore we can restrict - without loss of generality - to the class of 
centered spherical distributions, for which p. and S do coincide with the origin 
in M. k and the fc-dimensional identity matrix 1^, respectively. 

Under P&(0, 1/., /), the radial function f determines the distribution of ||X||. 
More precisely, ||X|| has probability density function fk(r) := (fJ,k-i;f)~ 1 i" k ~ 1 f(r) 
I[r>o] {I A stands for the indicator function of the Borcl set A); denote by Fk the 
corresponding distribution function. 

The assumption that fik-ij < oo guarantees that (2.1) is a density. The classi- 
cal Gaussian procedure for testing multivariate independence - namely, Wilks [34] 's 
test - requires the underlying distribution to have a finite variance (note that Puri 
and Sen [25] erroneously require finite moments of order four -an error which is 
repeated in Muirhead and Waternaux [20]). Consequently, when considering AREs 
with respect to Wilks' test, we will restrict to radial functions satisfying the stronger 
condition fJ.k+1-j '■= Jo°° r k+1 f(r) dr < oo, under which the distribution Pfc(0, Ik, f) 
has finite second-order moments. One can associate with each radial function / the 
radial function type of f defined as the class {/„, a > 0}, where / a (r) := /(ar), for 
all r > 0. By affine-invariance, one could restrict to parameter values of the form 
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(/i,E, /) = (0,Ifc,/ ao ) for which the variances of the associated elliptical distribu- 
tions arc equal to However, it will be convenient in the sequel to consider all 
possible radial functions, so that we will only fix (/i,S) = (0, Ifc.) and let / range 
over its radial function type. Some extremely mild smoothness conditions on /- 
which wc throughout will assume to be fulfilled are required to derive AREs. We 
refer to Taskinen, Kankainen and Oja [31] for details. 

The radial function / is said to be Gaussian if and only if / = <f) a for some 
a > 0, where <fi(r) := exp(— r 2 /2). Under Pfc(0, Ifc, </>), the probability density of 
||X|| is 4> k (r) := {2 ( - k ~ 2 ^ 2 T(k/2))~ 1 r k ~ 1 4)(r)I [r>0] where T(.) stands for the Euler 
gamma function; denote by the corresponding distribution function. Therefore, 

the distribution of ||X|j 2 = 1 (U)) 2 (throughout, U stands for a random variable 
uniformly distributed over (0,1)), still under Pfc(0, Ifc, </>), is xL an d the distribu- 
tion function of ||X|| is simply $fc(r) = ^(r 2 ), where denotes the chi-square 
distribution function with k degrees of freedom. 



2. 2. Testing for multivariate independence 

Consider an i.i.d. sample (X' 11; X 21 )', (X' 12 , X 22 )', ■ ■ ■ , (X' ln , X 2n )' of (p + ^-dimen- 
sional random vectors with the same distribution as (X^X^)'. We are interested 
in the null hypothesis under which Xi and X2, with clliptically symmetric distri- 
butions P p (fi 1 ,'Ei, f) and P q (u 2 , S 2 , <?), respectively, are mutually independent. As 
already mentioned, we can -without loss of generality (since the testing problem 
under study is invariant under block-diagonal affinc transformations) - restrict to 
centered spherical marginal distributions, that is, to assume (/i l5 Si) = (0,I p ) and 
(/x 2 ,E 2 ) = (0,I 9 ). 

The standard parametric procedure is Wilks [34] 's Gaussian likelihood ratio test 
{4>j\f, say), which rejects the null hypothesis (at asymptotic level a) whenever 



-nlog 



Sn S22 



> X 



y>q,l— a' 



where 



Sn S12 
S21 S22 



RVC; 



Xi, 

X2i 



X2i 



X 



stands for the partitioned sample covariance matrix, and where \ 2 



pq,l — a 



denotes the 



a-upper quantile of the chi-square distribution with pq degrees of freedom. 

The Taskinen, Kankainen and Oja [31] rank-based competitors of Wilks' test are 
defined as follows. Define the standardized subvectors Zn, . . . , Zi„ associated with 

~-l/2 

original subvectors Xn, . . . , Xi„ as Zn := Ti 1 (Xij — fa), i = 1, . . . , n, where fa 
and Ei are affinc-cquivariant root-n, consistent estimators of location and scatter, 
respectively. Consider the corresponding standardized spatial signs TJu = Zij/||Zij|| 
and let Ru denote the rank of ||Zij|| among ||Zn|j, . . . , ||Zi„|j; note that the scatter 
matrix estimate Si can be replaced by any affine-equivariant root-n consistent 
estimator of the shape matrix, as only the directions and ranks of distances are 
used in the analysis. The statistics U2i and i?2i are defined in the same way within 
the sample X 2 i,...,X 2 n- Letting K\^K 2 : (0,1) — > R be two square-integrable 
score functions, the (K\ , K2)-scoYe version of the rank test statistics for multivariate 
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independence proposed in Taskinen, Kankainen and Oja [31] is 

(2 . 2) ^ : = J ave . K2 Dl( „-,} 1 2 , 

where cr^- := J„ K 2 {u) du, and ||A|j 2 := tr(AA') is the squared Frobenius norm 
of A. Under the null hypothesis of independence (with elliptical marginals), this 
rank score statistic is asymptotically chi-square with pq degrees of freedom, and 
the associated test 4 , k 1 ,k 2 rejects the null hypothesis (still at asymptotic level a) 
as soon as T Ki ,k 2 > Xpq,i- a - 

As we show in the sequel, two particular cases (corresponding to two specific 
types of score functions) of the above rank score tests exhibit a remarkably good 
uniform efficiency behavior. 



3. Pitman non-admissibility of Wilks' test 



The van der Waerden (normal-score) version of the rank-based test statistic (2.2) is 
obtained with the score functions K\ = $~ 1 = 1 ) 1/2 , K 2 = <&„ 1 = ("J" 1 ) 1 / 2 : 



r, 



vdW 



i" 1 



i>-i 
i ' q 



Ro 



1 



u« u' 



In order to compute asymptotic relative efficiencies of the resulting van der Waer- 
den test 4> v dw with respect to Wilks' test we must embed the null hypothesis of 
independence in a model allowing for local (contiguous) alternatives of dependence. 
Such an embedding is not as obvious as in classical parametric models, where lo- 
cal values of the parameter usually provide the required alternatives. As in Gicscr 
and Randies [3] and Taskinen, Kankainen and Oja [31], we adopt (a multivariate 
extension of) the local alternatives considered by Konijn [19], of the form 



(3.1) 



(l-n-^Ip 



Y 2l 



1. 



where Yu and *Y 2 i denote mutually independent random vectors, with elliptic 
distributions P p (0, l p , f) and P ? (0, l q ,g), respectively, M is a non-random matrix 
with appropriate dimensions, and 5 G M. In that local model, the null hypothesis 
of independence reduces to Hq : 5 = 0. As shown in Taskinen, Kankainen and 
Oja [31], the asymptotic relative efficiency of the rank score test 4>k 1 ,k 2 with respect 
to Wilks' test does not depend on M. Their results imply that the asymptotic 
relative efficiency of the van der Waerden test 4> v dw based on T v dw with respect 
to Wilks' test, under the sequence of local alternatives in (3.1) is 



(3.2) KBJ£P q 'f g {<j) vd w/M 



1 



Ap 2 q 2 



D p {<t>, f)C q {cj>, g) + D q (<p, g)C p (0, /) 



where, denoting by <Pf{r) := — /' (r)/ f(r) the optimal location score function asso- 
ciated with the radial function /, we let 



Ckfaf) :=E *?{U)<p,{Ft\U)) and D^J) := E ^ \U) F^ l (U) 



190 



M. Hallin and D. Paindaveine 



Table 1 

AREs of the van der Waerden rank score test <j>vdw S or multivariate independence with respect 
to Wilks' test <f>j\f, under standard multivariate Student (with 3, 4, 6, and 12 degrees of 
freedom) and standard multinormal densities, for subvector dimensions p = 2 and 
q = 1,2, 3, 4, 6, and 10, respectively 



Up 


q 


Un 


3 


4 


6 


12 


oo 


1 


3 


1.378 


1.295 


1.266 


1.281 


1.339 




4 


1.293 


1.190 


1.141 


1.135 


1.167 




6 


1.267 


1.144 


1.078 


1.054 


1.067 




12 


1.285 


1.141 


1.058 


1.019 


1.016 




oc 


1.343 


1.174 


1.072 


1.017 


1.000 


2 


3 


1.400 


1.311 


1.277 


1.289 


1.343 




4 


1.311 


1.204 


1.152 


1.144 


1.174 




6 


1.277 


1.152 


1.085 


1.060 


1.072 




12 


1.289 


1.144 


1.060 


1.021 


1.017 




oo 


1.343 


1.174 


1.072 


1.017 


1.000 


3 


3 


1.417 


1.323 


1.286 


1.294 


1.346 




4 


1.325 


1.214 


1.161 


1.150 


1.179 




6 


1.286 


1.159 


1.091 


1.065 


1.076 




12 


1.292 


1.146 


1.062 


1.023 


1.019 




oo 


1.343 


1.174 


1.072 


1.017 


1.000 


1 


3 


1.430 


1.332 


1.292 


1.298 


1.348 




4 


1.336 


1.223 


1.167 


1.156 


1.183 




6 


1.294 


1.165 


1.096 


1.069 


1.080 




12 


1.295 


1.149 


1.064 


1.024 


1.020 




oo 


1.343 


1.174 


1.072 


1.017 


1.000 


6 


3 


1.448 


1.345 


1.301 


1.304 


1.351 




1 


1.353 


1.236 


1.177 


1.163 


1.189 




6 


1.306 


1.175 


1.103 


1.075 


1.085 




12 


1.300 


1.153 


1.068 


1.027 


1.023 




oo 


1.343 


1.174 


1.072 


1.017 


1.000 


10 


3 


1.471 


1.361 


1.312 


1.311 


1.353 




4 


1.375 


1.252 


1.190 


1.173 


1.196 




6 


1.323 


1.188 


1.114 


1.084 


1.092 




12 


1.308 


1.159 


1.073 


1.032 


1.027 




oo 


1.343 


1.174 


1.072 


1.017 


1.000 



Some numerical values of these AREs, under multivariate t- and normal distribu- 
tions, are provided in Table 1. All these values are larger than or equal to 1, and 
seem to be equal to 1 only if both marginals are multinormal. 

Taskinen, Kankainen and Oja [31] point out that it is remarkable that, in the 
multinormal case, the limiting Pitman efficiency of the van der Waerden score test 
4>vdw equals that of Wilks' test. But, it is even more remarkable that, as we shall 
see, the multinormal case is actually the least favorable one for the van der Waerden 
procedure. Proposition 1 below indeed states that, as soon as one of the marginals is 
not Gaussian, the van der Waerden test strictly beats Wilks' test (Table 1 provides 
an empirical confirmation of this fact). The Pitman non-admissibility of Wilks' test 
follows. 

Proposition 1. For all integers p, q > 1 and all radial functions f,g such that 
Mp+i;/ < 00 an d fJ-q+i-.g < °°j we have 

ARW q f g (<j> vdw /M > 1, 

where equality holds iff f and g are Gaussian, with common scale. 

To prove this proposition, we need the following intermediate result; see Pain- 
daveine [22] for an elementary proof, based on an idea developed by Gastwirth and 
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Wolff [2]. 

Lemma 1. For all integer k > 1 and all radial function f satisfying fJ<k+i;f < oo, 
we have Dk(4>, f)C k (4>, f) > k 2 , where equality holds iff f is Gaussian. 

Proof of Proposition 1. The proof is based on the decomposition 

(D p {4>, f)C q (<P, g) + D q (d>, g)C p {4>, /)) * = A v q J g + B**, 

where 

A*f := 4 D p (<f>, f)C p {<j>, f)D g (<f>, g)C q {<j>, g) 

and 

B*t := (D p {cj>J)C q {cj>,g)-D q {<j>,g)C p {cj>J)) 2 . 
It directly follows from Lemma 1 that A^J > 4p 2 q 2 , so that (3.2) entails 

(3.3) ARE^(W<^) = ^(^i' + B»t) > Atf > 1. 

Let us now show that equality holds iff / and g are Gaussian, with the same vari- 
ance. Equality in (3.3) indeed requires A?'? = 4p 2 q 2 and B^'f = 0. From Lemma 1, 
A v q ^ g = Ap 2 q 2 implies that both / and g are Gaussian (/ = <f> a and g = say). Now, 
since D k ((f),(j) a ) = or 1 D k {4>,4>) = a~ 1 k and C k ((f),(j) a ) = aC k {<j),cj>) = aD k {(f>,4>) = 
ak for all k, we have B*'%° = p 2 q 2 ((b/a) — (a/b)) 2 , which is equal to zero iff a = b. 
Consequently, equality holds iff / = g = <p a , for some a > 0. □ 



4. A Hodges-Lehmann result for multivariate independence 

We now turn to the Wilcoxon test statistic 

Tw ~ II ave * i R u U « U 2J II 2 > 

which is associated with the score functions Ki(u) = K%{u) = u, u <G (0, 1). The 
asymptotic relative efficiency of the corresponding Wilcoxon test <pw with respect 
to Wilks' test </>m, under the sequence of local alternatives (3.1), is given by 

ARE^^w/M = ^-(D p {IJ)C q {I,g) + D q (I,q)C p {IJ))\ 

where we let 

C k (I,f) ^^[UipfiF^iU))] and D k (I,f) := E[uFj^\U) 

see Taskinen, Kankainen and Oja [31]. Some numerical values of these AREs, under 
multivariate t- and normal distributions, are provided in Table 2. The uniformly 
good asymptotic efficiency of the Wilcoxon test in Table 2 is not just an empirical 
finding, as shown by the following result which provides the lower bound of these 
AREs for any fixed values of the dimensions p, q of the marginals (some numerical 
values of this lower bound are presented in Table 3). 
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Table 2 

AREs of the Wilcoxon rank score test <fiw for multivariate independence with respect to Wilks' 
test <f>j*f, under standard multivariate Student (with 3, 4, 6, and 12 degrees of freedom) 
and standard multinormal densities, for subvector dimensions p = 2 and q = 1, 2, 3, 4, 6 

and 10, respectively 



Up 


q 


Un 


3 


4 




6 


12 


oo 


1 


3 


1.305 


1.227 




1.193 


1.193 


1.222 




4 


1.239 


1.147 




1.099 


1.085 


1.098 




6 


1.208 


1.104 




1.044 


1.018 


1.018 




12 


1.204 


1.086 




1.015 


0.978 


0.969 




oo 


1.219 


1.087 




1.006 


0.959 


0.940 


2 


3 


1.305 


1.237 




1.211 


1.219 


1.257 




1 


1.237 


1.152 




1.111 


1.102 


1.121 




6 


1.211 


1.111 




1.055 


1.033 


1.037 




12 


1.219 


1.102 




1.033 


0.997 


0.989 




oo 


1.257 


1.121 




1.037 


0.989 


0.970 


3 


3 


1.274 


1.213 




1.193 


1.206 


1.248 




4 


1.203 


1.125 




1.089 


1.084 


1.106 




6 


1.179 


1.084 




1.032 


1.013 


1.020 




12 


1.192 


1.079 




1.013 


0.980 


0.973 




oo 


1.245 


1.110 




1.027 


0.979 


0.960 


1 


3 


1.248 


1.192 




1.175 


1.191 


1.235 




4 


1.174 


1.101 




1.068 


1.066 


1.090 




5 


1.149 


1.059 




1.011 


0.993 


1.002 




12 


1.165 


1.056 




0.992 


0.961 


0.955 




oo 


1.228 


1.095 




1.013 


0.966 


0.947 


6 


3 


1.211 


1.161 




1.150 


1.168 


1.215 




1 


1.134 


1.067 




1.038 


1.039 


1.066 




6 


1.105 


1.022 




0.978 


0.963 


0.974 




12 


1.122 


1.019 




0.959 


0.930 


0.927 




oo 


1.198 


1.068 




0.988 


0.943 


0.924 


10 


3 


1.173 


1.129 




1.121 


1.144 


1.193 




4 


1.090 


1.029 




1.005 


1.009 


1.038 




6 


1.056 


0.979 




0.940 


0.929 


0.941 




12 


1.069 


0.973 




0.918 


0.892 


0.891 




oo 


1.158 


1.033 




0.955 


0.911 


0.893 








Table 3 








Some numerical values, for various values of the dimensions p and q 


of the subvectors, 


of the 


Hodges-Lehmann lower bound for the asymptotic relative efficiency 


of the Wilcoxon rank 


s 


core test 


4>w f or multivariate independ 


ence with respect to Wilks ' test <f>pj- 




p/q 


1 


2 


3 


4 


6 


10 


oo 


1 


0.856 


0.884 


0.867 


0.850 


0.826 


0.797 


0.694 


2 




0.913 


0.895 


0.878 


0.853 


0.823 


0.717 


3 






0.878 


0.861 


0.836 


0.807 


0.703 


4 








0.845 


0.820 


0.792 


0.689 


6 










0.797 


0.769 


0.669 


10 












0.742 


0.646 


oo 














0.563 


Proposit 


ion 2. 


Let p, q > 1 be two integers. Th 


en, letting 








Cfe : 


= inf |x > 




1/2 ( x 


)'-»}• 


k £No, 





where J r denotes the first-kind Bessel function of order r, the lower bound for 
the asymptotic relative efficiency of <pw with respect to <f>jy, for fixed subvector 
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dimensions p, q, is 

MARE^w/M = (2cl +p if (2 4 + q if , 

where the infimum is taken over the collection of radial functions /, g for which 
fJ-p+i-J < o° an d fJ-q+i-.g < 00. The infimum is reached at 

(f>g)€ {(h P A r )> h <iA r )) : = (h p ,i{ar), h q ,i (arn , a > oj , 

where hk,i denotes "the" (7ife,i is defined up to a positive scalar factor) radial func- 
tion associated with the radial cumulative distribution function 

„ l \ Vr J V2k=T/2 ( r ) T T 

H k,l(r) := =-= j— t i[o <r <c k ] + J [r>c fc ]- 

To prove this proposition, we need the following lemma, which is established in 
the proof of Proposition 7 of Hallin and Paindaveinc [is] . 

Lemma 2. Let k > 1 &e a /ixec? integer. Then, 



M{ Dk (IJ)C k (IJ)} = ^- 



where the infimum is taken over the collection of radial functions f for which 
Mfc+i;/ < °°; an ^ ^ e infimum is reached at the radial functions f £ {/ifc jCr (r), cr > 
0}. Moreover, letting u>k '■= (2c^ + fc — l)/(8c/-) } we /iai>e Dk{I,hk )U>k ) = 1. 

Since we have D k (I,f a ) = a~ 1 D k (I,f) and C k (I,f a ) = aC k (I,f) for all fc, 
the quantity D k (I , f a )C k (I , f a ) does not depend on a. This allows for identifying 
a particular member /ife i<T( . of the radial function type {h kt(7 (r), a > 0} such that 
D k (I , h kMk ) = 1. According to Lemma 2, a k = cu k . 

Proof of Proposition 2. Proceeding as in the proof of Proposition 1, we consider 
the decomposition 

where 

A ll ~-D p (I,f)C p (I,f)D q (I,g)C q (I,g), 



pq 

B li ■■= J- (D p (I,f)C q (I,g)-D q (I,g)C p (IJ)) 



4pq 

Lemma 2 directly yields that, for all couple (f,g) of radial functions, 

(4.1) ARE^/M > AH > zU^T2 ( 2c p +P~ !) 2 ■ 

We now show that the right-hand side in (4.1) actually coincides with the in- 
fimum, by determining the (non-empty) collection of couples (/, g) achieving the 
bound in (4.1). At the couple (/, g), the bound is achieved iff 

( 4 - 2 ) AP d = 2 Toi ? ( 2 4+P- 1 ) 2 ( 2 ^ + ^ 1 ) 2 I and 

* pq^p^q 

B pJ = 
^qg u - 
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Lemma 2 shows that (4.2) holds iff 

ft g> i(&r)) J a,&>o}. 
Now, using the fact that Dk(I, hk,i) = for all k, we have 

KtZ = 4^ ( \ D ?( J > VOQ^: Kx) - \ D q (I, h q!l )C p (I, Vi)) 2 

= A ( ^D q (I, h q<1 )C q (I, h g ,i) - ^D P (I, h Ptl )C p (I, h p ,!)) 2 
= 

iff 

(UqO\ 2 = D q (I 7 h qA )C q (I 7 h qA ) = c 2 p (2c 2 q + g~l) 2 = 
Vw^J D p (I,h pA )C p (I,h pA ) c 2 {2c 2 + P ~l) 2 \ij p ) ' 
that is, iff a = b. □ 

Somewhat surprisingly, the main results of this paper are scale-dependent. Es- 
sentially, the bounds in Propositions 1 and 2 are achieved at densities / and g with 
common scale. This at first sight is puzzling. Indeed, all tests involved are (block-) 
affmc-invariant. This dependence on scale however is entirely due to the form of 
the local Konijn alternatives (3.1) considered throughout. Appropriate rescaling 
of M would take care of it. For the sake of coherence with the literature (see, 
e.g., Giescr and Randies [•">], Taskincn, Kankaincn and Oja [30, 3f] and Taskincn, 
Oja and Randies [32]), however, we sticked to the traditional definition of Konijn 
alternatives. 

5. Final comments 

Multivariate signed rank tests, based on affine-invariant concepts of signs and 
ranks -Randies' interdirections (Randies [28]); (pseudo)-Mahalanobis ranks (Hallin 
and Paindaveine [7, 8, 9, fO, 11]; Hallin, Oja and Paindaveine [6]); Tyler's angles 
(Hallin and Paindaveine [9]); Oja-Paindaveine ranks (Oja and Paindaveine [21])- 
recently have attracted renewed attention to rank-based methods in the context of 
multivariate analysis, a domain where the 1971 monograph by Puri and Sen has 
remained the main reference for more than thirty years. 

The specific problem of testing independence between elliptical random vectors 
has been investigated in a series of papers by Taskincn, Kankaincn and Oja [30, 31] 
and Taskincn, Oja and Randies [32]. Simulations show that the performances of 
the procedures proposed in these papers are quite good. This empirical finding 
is confirmed here by an investigation of their asymptotic relative efficiencies with 
respect to their traditional Gaussian counterpart, Winks' test of independence, based 
on classical correlations. More specifically, we obtain, for the Wilcoxon and van der 
Waerden versions of the Taskinen-Kankainen-Oja tests, analogues of the classical 
Chcrnoff- Savage [1] and Hodges-Lehmann [15] results. The Chcrnoff-Savage result 
of Proposition 1 in particular establishes the non-admissibility of Wilks' procedure, 
which is uniformly dominated, in the Pitman sense, by the van der Waerden version 
of Taskincn, Kankaincn and Oja [31]. 
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